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We present measurements of net charge fluctuations in Au + Au coflisions at ^^snn = 19.6, 62.4, 
130, and 200 GeV, Cu + Cu collisions at ^/s]^ = 62 A, 200 GeV, and p + p collisions at = 200 
GeV using the dynamical net charge fluctuations measure — ,dyn. We observe that the dynamical 
fluctuations are non-zero at all energies and exhibit a modest dependence on beam energy. A weak 
system size dependence is also observed. We examine the collision centrality dependence of the 
net charge fluctuations and find that dynamical net charge fluctuations violate 1/Nch scaling, but 
display approximate 1/Npart scaling. We also study the azimuthal and rapidity dependence of the 
net charge correlation strength and observe strong dependence on the azimuthal angular range and 
pseudorapidity widths integrated to measure the correlation. 

PACS numbers: 25.75.Gz, 25.75.Ld, 24.60. Ky, 24.60.-k 

Keywords: Net charge fluctuations, azimuthal correlations, QGP, Heavy Ion Collisions 



INTRODUCTION cators of the formation of a quark gluon plasma (QGP) in 



Anomalous transverse momentum and net charge 
event-by-event fluctuations have been proposed as indi- 
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the midst of high-energy heavy ion coUisions. A number 
of authors P, U, 0] have argued that entropy conserving 
hadronization of a plasma of quarks and gluons should 
produce a final state characterized by a dramatic reduc- 
tion of the net charge fluctuations relative to that of a 
hadron gas. Simply put, their prediction relics on the 
notion that quark-quark correlations can be neglected, 
and hadronization of gluons produces pairs of positive 
and negative particles not contributing to the net charge 
fluctuations. Accounting for the fractional charge of the 
quarks, they find that, for a QGP, the variance of the 
ratio of positive and negative particles scaled by the to- 
tal charged particle multiplicity, a quantity they call D, 
should be approximately four times smaller than for a 
gas of hadron. Quark-quark correlation may however 
not be negligible; Koch et al. [l| extended their origi- 
nal estimates to include susceptibilities calculated on the 
lattice. They find the quantity D = 4:{AQ^/Nch (where, 
AQ^ is the variance of the net charge, Q = - N~ 
and Nch is the total number of charged particles observed 
in a particular momentum space window under consid- 
eration) is quantitatively different from their first basic 
estimate but nonetheless still dramatically smaller than 
values expected for a hadron gas. It is clear however 
that hadron collisions, and in particular heavy ion col- 
lisions, produce substantial number of many high mass 
particles, and specifically short lived (neutral) particles 
or resonances which decay into pairs of positive and neg- 
ative particles. Such decays increase the multiplicity of 
charged particles in the final state while producing negli- 
gible impact on the net charge variance. Jeon and Koch 
have in fact argued that one can use the magnitude of net 
charge fluctuations to estimate the relative production 
of p and uj mesons Q . Calculations based on a thermal 
model lead to a value D of order of 2.8 which although 
reduced relative to the value expected for a pion gas is 
nonetheless remarkably larger than that predicted for a 
QGP Note that transport models such as UrQMD 
predict values in qualitative agreement with those of ther- 
mal models [sl] . A measurement of net charge fluctuations 
therefore appears, on the outset, as an interesting means 
to identify the formation of quark gluon plasma in high- 
energy heavy ion collisions. 

First measurements of net charge fluctuations were re- 
ported by both PHENIX 6] and STAR J] collaborations 
on the basis of Au + Au data acquired during the flrst 
RHIC run at ^/sNN = 130 GeV. Measurements were re- 
ported by PHENIX ^ in terms of a reduced variance, 
ujQ = {AQ'^ / Nch- Unfortunately, measured values of 
this quantity depend on the efficiency. STAR instead 
reported results [3, [1| in terms of dynamical net charge 
fiuctuations measure, _dyn, which is found to be a ro- 
bust observable i.e., independent of detection efficiency, 
z^-i .dyn is defined by the expression: 

_ {N+iN+ - 1)) , (A^_(jV_ - 1)) ^ {N^N+) 

(1) 



where N± are the number of positively and negatively 
charged particles in the acceptance of interest. Note 
that there exists a simple relationship between uq and 
'^+-Ayn written J^+-,dyri = 4(wQ - l)/Nch- This rela- 
tionship is however applicable only if ujq is corrected for 
finite detection effects. Because such corrections are not 
trivial, we favor the use of i>-\ .dyn- We note addition- 
ally that both UJQ (corrected for efficiency) and vjf ^dyn 

may be expressed (at least approximately) in terms of 
the variable D ~ Nch{AR^ (with R = N+/N^) used 
by Koch et al. [ij for their various predictions. Their 
use, for experimental measurements, avoid pitfalls asso- 
ciated with measurements of average values of the ratio, 
R, of particle multiplicities — where the denominator, 
N- may be small or even zero The measurements 
performed by the STAR 7] and PHENIX [6] collabora- 
tions showed the dynamical net charge fluctuations in 
Au + Au at y^SNN = 130 GeV are finite but small rel- 
ative to the predictions by Koch et al. [l| for a QGP. 
The magnitude of the net charge fluctuations was found 
to be in qualitative agreement with HIJING predictions 
although the data exhibit centrality dependence not 
reproduced within the HIJING calculations. Measured 
values also qualitatively agree with predictions by Bialas 
for quark coalescence [ifl] and Koch et al. for a resonance 
gas [1]. 

The scenario for dramatically reduced fluctuations as 
an evidence for the formation of QGP is clearly excluded 
by the data at 130 GeV. However, in light of predictions 
of a tri-critical point of the equation of state in the range 
10 < y/sNN < 60 GeV [m, [13], one might argue the re- 
duction of fluctuation might be larger at lower beam ener- 
gies. Conversely, one may also argue the volume of QGP 
formed in Au + Au collisions might increase at higher 
beam energies leading to reduced fluctuations at higher 
beam energy instead. One is thus led to wonder whether 
the fluctuations may be found to vary with beam energy 
thereby indicating the production of QGP above a criti- 
cal threshold, or with progressively increasing probability 
at higher energies. In this paper, we consider this possi- 
bility by investigating how the strength of the dynamical 
net charge fluctuations vary with beam energy and sys- 
tem size in Au + Au and Cu + Cu collisions ranging in 
center of mass energy from the highest SPS energy to 
the highest RHIC energy, and relative to p + p collisions 
at y^SNN = 200 GeV. The analysis presented also pro- 
vides, independent of existing models, new information 
that may shed light on the collision dynamics. 

Various issues however complicate the measurement 
and interpretation of net charge fluctuations. First, one 
must acknowledge that particle flnal state systems pro- 
duced in heavy ion collisions although large, are nonethe- 
less flnite and therefore subject to charge conservation 
effects. Produced particles are also measured in a flnite 
detector acceptance. Second, one may question whether 
the dynamical net charge fluctuations produced within 
the QGP phase may survive the hadronization process 
[l^. Shuryak and Stephanov have argued based on 
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solutions of the diffusion equation within the context of a 
model involving Bjorken boost invariance, that diffusion 
in rapidity space considerably increases the net charge 
fluctuations. They further argued that the reduced fluc- 
tuations predicted for a QGP might be observable only if 
fluctuations are measured over a very large rapidity range 
(of order of 4 units of rapidity). Unfortunately, charge 
conservation effects increase with the rapidity range con- 
sidered and might become dominant for rapidity ranges 
of four units or more. Gavin et al. however argued 
that the classical diffusion equation yields non-physical 
solutions in the context of relativistic heavy ion collisions. 
They proposed a causal diffusion equation as a substi- 
tute of the classical diffusion equation for studies of net 
charge fluctuation dissipation. They found that causal- 
ity substantially limits the extent to which diffusion can 
dissipate these fluctuations. 

Third, there exists the possibility that the treatment 
by Koch et al. [l| of quark and gluons behaving as inde- 
pendent particles carrying full entropy may be inappro- 
priate. Consider for instance that recent measurements 
of elliptical anisotropy of particle emission in Au -\- Au 
collisions show that meson and baryon elliptical flow, W2, 
scales in proportion to the number of constituent quarks 
for transverse momenta in the range 1-4 GeV/c, thereby 
suggesting that hadrons are produced relatively early in 
the collisions through "coalescence" or recombination of 
constituent quarks. In a constituent quark scenario, the 
role of gluons in particle production is reduced. Rela- 
tively smaller charged particle multiplicities are therefore 
expected, and net charge fluctuations are correspond- 
ingly larger. Bialas [1^ conducted a simple estimate of 
such a scenario, and reported net charge fluctuations D 
may be of order 3.3. Interestingly, this estimate suggests 
fluctuations might be even larger than that expected for 
a resonance gas, and as such should also be identifiable 
experimentally. 

Theoretical estimates of the effect of hadronization on 
net charge fluctuation have been for the most part re- 
stricted to studies of the role of resonances, diffusion 
[il m, [ll, and thermalization [13, d. One must 
however confront the notion that collective motion of 
produced particles is clearly demonstrated in relativistic 
heavy ion collisions. Voloshin pointed out [l^ that in- 
duced radial flow of particles produced in parton-parton 
collisions at finite radii in nucleus- nucleus collisions gen- 
erate momentum-position correlations not present in ele- 
mentary proton-proton collisions. Specifically, the effect 
of radial flow is to induce azimuthal correlations and to 
modify particle correlation strengths in the longitudinal 
direction. Voloshin showed that two-particle momentum 
correlations (ApTi,ApT2) are in fact sensitive to radial 
velocity profile as well as the average flow velocity. While 
one may not intuitively expect net charge fluctuations to 
exhibit a dramatic dependence on radial flow, simulations 
based on a simple multinomial particle production model 
including resonances such as the p(770) indicate that net 
charge correlations are in fact also sensitive to radial flow 



through azimuthal net charge correlations [20|. They 
may as such be used to complement estimates of radial 
velocity obtained from fits of single particle spectra with 
blast-wave parameterization or similar phenomenologies. 

Measurements of charged particle fluctuations have 
also been proposed as a tool to discriminate between 
predictions of various microscopic models of nuclear col- 
lisions. Zhang et al. [21] find that measurements of dy- 
namical fluctuations should exhibit sensitivity to rescat- 
tering effects based on calculations without rescattering 
with models VNIb ^ and RQMD They also found 
that models VNIb, HIJING [9], HIJING/BB [H and 
RQMD predict qualitatively different dependencies on 
collision centrality. Similar conclusions are obtained by 
Abdel-Aziz [25]. 

Bopp and Ranft [1^ compared predictions of net 
charge fluctuations (at mid rapidities) by the dual par- 
ton model and statistical (thermal) models, and found 
significant differences in the dispersion of the charges pre- 
dicted by these models. They hence argued that charged 
particle fluctuations should provide a clear signal of the 
dynamics of heavy ion processes, and enable a direct mea- 
surement of the degree of thermalization reached in heavy 
ion collisions. Gavi n [TtI [Tsj similarly argued, based on 
data by PHENIX [i,|23| and STAR 0,113 that measured 
transverse momentum and net charge fluctuations indeed 
present evidence for thermalization at RHIC. 

In this work, we present measurements of dynamical 
net charge fluctuations in Au -\- Au collisions at ^/snn = 
19.6, 62.4, 130, and 200 GeV, Cu + Cu collisions at 
y/sMN = 62.4, 200 GeV and inp+p collisions at -/s = 200 
GeV. We study the beam energy, system size and collision 
centrality dependencies quantitatively in order to iden- 
tify possible signature of the formation of a QGP. Some 
of the results presented in this work have been reported 
as preliminary data at various conferences [20| . The pa- 
per is organized into sections on Experimental Method, 
Results, Systematic Error Studies, and Conclusions. 



II. EXPERIMENTAL METHOD 

Our study of dynamical net charge fluctuations depen- 
dence on the beam energy is based on the observable 

(jyn used in the first STAR measurement The 

definition of _dyn, its properties, and relationships to 

other measures of event- by- event net charge fluctuations 
were motivated and presented in detail in Refs. (sl. |29|. 
The robustness of .dyn as an experimental observ- 
able was also discussed on the basis of Monte Carlo toy 
models by Nystrand et al. [sO]; the authors verified ex- 
plicitly with simple Monte Carlo generators that ^dyn 

is insensitive to the details of the detector response and 

efficiency. Indeed, they verified that values of dyn are 

independent of the track detection efficiency when the ef- 
ficiency is uniform over the measured kinematic range. If 
the efficiency is not perfectly uniform across the accep- 
tance, the robustness of dyn is reduced in principle. 
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However, in this work, the acceptance of the measure- 
ment is hmited to a kinematic range where the efficiency 
is essentially uniform, and such effects are, therefore, neg- 
ligible. 

We here briefly review the definition and essential 
properties of this observable. Rather than measuring the 
event-by-event fiuctuations of the ratio of positive and 
negative particle multiplicities (in a given acceptance), 
one considers the second moment of the difference be- 
tween the relative multiplicities (N^) and N^/ (-/V_) 
as follows 



/ N+ iV_ 



(2) 



The Poisson limit, ^stat of this quantity is equal to 



1 



l^A — ,stat 



{N+) (N.) 



(3) 



The "non-statistical" or "dynamical" fluctuations can 
thus be expressed as the difference between the above 
two quantities: 



— ,dyn 



- v+ 
{N+iN+ 



- ,stat 

-i)> 



(4) 



(jV_(jV_-l)) 



g {N+N_) 
(iV_) {N+) 

From a theoretical standpoint, ^dyn can be ex- 
pressed in terms of two-particle integral correlation func- 
tions as ^dyn = R++ + R — 2R^ , where the terms 

Rap are ratios of integrals of two and single particle pseu- 
dorapidity density functions defined as follows : 



Raf} 



J driadr]f3 



dN 



dN 
' drjct 



(5) 



The dynamical net charge fluctuations variable _dyn 

is thus basically a measure of the relative correlation 

strength of , and H — particles pairs. Note that 

by construction, these correlations are identically zero for 
Poissonian, or independent particle production. In prac- 
tice, however, produced particles are partly correlated, 
either through the production of resonances, string frag- 
mentation, jet fragmentation, or other mechanisms. The 

relative and absolute strengths of i?--, and R-i 

may vary with colliding systems, and beam energy. In ad- 
dition, by virtue of charge conservation, the production 
of H — pairs is expected to be more strongly correlated 

than the production of or pairs. For this reason, 

it is reasonable to expect to be larger than or 

R In fact, one finds experimentally that 2R^ is ac- 
tually larger than the sum + R in p + p and p + p 



collisions measured at the ISR and FNAL [3i|,l33. Mea- 
surements of dyn are thus expected and have indeed 

been found to yield negative values in nucleus- nucleus 
collisions as well [tJ. 

We also note i^^ ^dyn is essentially a measure of the 

variance of — iV_ . This difference is "orthogonal" to 
the multiplicity iV_|_-|-7V_, and thus linearly independent. 

There is, therefore, no bias introduced by binning ly^ ^dyn 

measurements on the basis of the reference multiplicity 
(multiphcity within \r]\ <0.5) as discussed below. 

As a technical consideration, our study of the i^^ dyn 

dependence on collision centrality is carried out in terms 
of charged particle multiplicity bins, as discussed in de- 
tail below. To avoid dependencies on the width of the 
bins, we flrst determine the values of dynamical fluc- 
tuation, ,dyn(w), for each value of multiplicity, m. 

The dynamical fluctuations are then averaged across the 
selected finite width of the centrality bins with weights 
corresponding to the relative cross section, p{m), mea- 
sured at each value of multiplicity. For example, in the 
multiplicity range from rrimin to nimax, we calculate the 
average as follows: 



— ,dyn (^min — ^ ^ ^max') 



Y,v+^^dyn{m)p{m) 



(6) 

This study is based on the notion that if Au + Au 
collisions (or any other A + A system) trivially consist 
of a superposition of independent nucleon-nuclcon colli- 
sions, with no rescattering of the produced secondaries, 
then _dyn is expected to scale inversely to the num- 
ber of participating nucleons and the number of charged 
particles, or more appropriately, the number of actual 
nucleon + nucleon collisions. One can thus infer that the 
quantity ,dyixdN ch / drj\ should be independent of col- 
lision centrality under such a scenario. We shall therefore 
examine whether indeed the dynamical net charge fiuc- 
tuations scale with the number of participants, or the 
invariant multiplicity. 

The data used in this analysis were measured using 
the Solenoidal Tracker at RHIC (STAR) detector dur- 
ing the 2001, 2002, 2004 and 2005 data RHIC runs at 
Brookhaven National Laboratory. They include Au ^ Au 
coUisions data collected at ^/snn = 19.6, 62.4, 130, and 
200 GeV, Cu + Cu collisions data at ^/smn — 62.4, 
200 GeV and p + p collisions data measured at i/i = 
200 GeV. The Au + Au coUisions at y^sWJ^ = 62.4, 130, 
and 200 GeV data and Cu+Cu collisions at ^/snn = 62.4 
and 200 GeV were acquired with minimum bias triggers 
accomplished by requiring a coincidence of two Zero De- 
gree Calorimeters (ZDCs) located at 18 m from the center 
of the interaction region on either side of the STAR de- 
tector. For 19.6 GeV data, a combination of minimum 
bias and central triggers was used. The centrality trigger 
was achieved using a set of scintillation detectors, called 
the Central Trigger Barrel (CTB) surrounding the main 
Time Projection Chamber (TPC). Technical descriptions 
of the STAR detector and its components are published 
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in technical reports [s^ Is^l ■ For p + p collisions a mini- 
mum bias trigger was used based on the CTB detector. 
The analysis carried out in this work is rather similar to 
that published in the first net charge fluctuation mea- 
surement 7\. 

This analysis is based on charged particle track recon- 
struction measurements performed with the STAR-TPC. 
The TPC is located in a large solenoidal magnetic field 
producing a uniform axial magnetic field. The magnetic 
field was set to 0.25 T for Au + Au collisions at y^sjvw = 
19.6 and 130 GeV data, and 0.5 T for Au + Au and 
Cu -t- Cu collisions at 62.4, and 200 GeV data. The in- 
creased magnetic field results in a slight reduction of the 
detection efficiency for charged particle tracks with trans- 
verse momenta below 0.2 GeV/c, and a modest improve- 
ment in momentum resolution. This analysis used tracks 
from the TPC with transverse momentum in the range 
0.2 < pT < 5.0 GeV/c with pseudorapidity I77I < 0.5. 
Systematic effects associated with finite thresholds and 
momentum dependent efficiency are discussed in Section 
IV. 

In order to limit the net charge fluctuations analysis 
to primary charged particle tracks only (i.e. particles 
produced by the collision), tracks were selected on the 
basis of their distance of closest approach (DCA) to the 
collision vertex. DCA is defined as the distance between 
the track and the primary vertex position. A nominal 
cut of DCA < 3 cm was used for results presented in 
this paper. Systematic effects associated with this cut 
are discussed in Section IV. 

Events were selected for analysis if their collision ver- 
tex lay within a maximum distance from the center of the 
TPC and they passed a minimum track multiplicity cut 
(see below). The vertex position was determined using 
a fit involving all found tracks. The maximum distance 
along the beam axis from the center of the TPC (also 
called the z vertex cut) was set to 75 cm for the Au + Au 
19.6 and 130 GeV data, further restricted to 25 cm for 
62.4 and 200 GeV data. However, a z vertex cut of 30 cm 
was used in Cu+Cu 62 A and 200 GeV data. A maximum 
of 75 cm was used for the p + p data. The wide 75 cm 
cut was used to maximize the event sample used in this 
analysis. The observable ^dyn measured in this analy- 
sis (as described below) is a robust experimental variable, 
and is by construction largely insensitive to restricted de- 
tection efficiencies provided those efficiencies do not vary 
dramatically across the detector acceptance. We indeed 
find that as long as the longitudinal cut is limited to val- 
ues below 75 cm, for which the track detection efficiency 
is rather insensitive to the pseudorapidity of the track, 

the measured values of dyn sue invariant within the 

statistical uncertainties of the p+p measurements. With 
a larger cut, the efficiencies drop dramatically at large 
rapidities, and ^dyn exhibits somewhat larger devia- 
tions. The analyses reported in this paper are based on 
100k, IM, 144k, lOM Au + Au events at 19.6, 62, 130, 
and 200 GeV, respectively, 9M and 5.5M Cu + Cu events 
at 62 and 200 GeV, and 2.7M p -fp events. 



The magnitude of net charge fluctuations is quite ob- 
viously subject to change with the total multiplicity of 
produced charged particles. It is thus necessary to mea- 
sure the magnitude of the fluctuations and correlations 
as a function of the collision ccntrality. Measurements 
at the AGS, SPS, and RHIC have shown that there is 
a strong anti correlation between the number of collision 
spectators (i.e. projectile/target nucleons undergoing lit- 
tle or no interaction with target/projectile nucleons) and 
the multiplicity of charged particles produced in the col- 
lisions. We use the standard collision ccntrality defini- 
tion used in other STAR analyses and base estimates 
of the collision ccntrality on the uncorrected multiplic- 
ity of charged particle tracks measured within the TPC 
in the pseudorapidity range -0.5 < rj < 0.5. While low 
multiplicity events correspond to peripheral (large im- 
pact parameter) collisions, high multiplicities are asso- 
ciated with central (small impact parameter) collisions. 
The pseudorapidity range -0.5 < rj < 0.5 is used for col- 
lision ccntrality estimates rather than the full range - 
1.0 < 77 < 1.0 in principle measurable with the TPC, 
to minimize effects of detector acceptance and efficiency 
on the collision ccntrality determination. With the nar- 
row cut - 0.5 < T] < 0.5, the track detection efficiency is 
rather insensitive to the position of the collision vertex 
(along the beam direction) in the range — 75 < 2; < 75 
cm used in the analysis of Au + Au at ^/snn = 130 GeV 
data, and ccntrality selection biases are thus negligible. 
The efficiency for tracks with 0.5 < < 1 on the other 
hand drops markedly for vertex positions \z\ > 50 cm. 
Although the analysis of Au + Au and Cu + Cu collisions 
at ^snn = 62.4 and 200 GeV data were conducted with 
the narrower \z\ < 25 cm and \z\ < 30 cm range, respec- 
tively, enabled by the more compact interaction region 
delivered by the accelerator during these runs, the ccn- 
trality determination was estimated on the basis of the 
same pseudorapidity range in order to provide uniform 
and consistent centrality cuts. 

The centrality bins were calculated as a fraction of this 
multiplicity distribution starting at the highest multiplic- 
ities. The ranges used were 0-5% (most central colli- 
sions), 5-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 
60-70%, and 70-80% (most peripheral) for Au + Au colli- 
sions. Similarly, collision centrality slices used in Cu+Cu 
collisions are 0-10% (most central), 10-20%, 20-30%, 30- 
40%, 40-50% and 50-60% (most peripheral). Each cen- 
trality bin is associated with an average number of partic- 
ipating nucleons, Npart, using Glauber Monte Carlo cal- 
culation [sH]. At low multiplicities, the finite detector ac- 
ceptance and track detection efficiencies imply estimates 
of the collision centrality are subject to large errors. 

Events arc included or "counted" in this analysis pro- 
vided a collision vertex is found (as per the discussion 
of the previous paragraphs) and at least one particle is 
found in the range —0.5 < rj < 0.5. While event count- 
ing efficiencies are essentially unity for large multiplicity 
collisions, they are limited (< 1) for small multiplicities 
corresponding to most peripheral collisions. The limited 
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efficiency stems from finite track and vertex finding ef- 
ficiencies. Track finding efficiency within the TPC was 
studied through detailed Monte Carlo simulations of the 
detector response with track embedding. For minimal 
track quality cuts such as those used in this analysis, 
one finds the track finding efficiency is of order 95% for 
Pt > 0.2 GeV/c in peripheral collisions. It reduces to 
approximately 85% for most central collisions and falls 
to zero for primary tracks with px < 0.1 GeV/c. The 
efficiencies of positive and negative particles are found to 
be the same within the statistical errors. The data shown 
were integrated for tracks with 0.2 < px < 5.0 GeV/c, 
\r]\ < 0.5 and < (j) < 27r. Note that the minimum cut 
used in this new analysis is different than that used in the 
first reported study 7]. A value of 0.2 GeV/c is used for 
all measured beam energies and field settings to avoid 
systematic effects associated with px dependent detec- 
tion efficiency below 0.2 GeV/c. The results presented in 
this work for 130 GeV are nonetheless in agreement with 
results reported by STAR in the first measurement of net 
charge^fluctuations in Au + Au collisions at ^/snn = 130 
GeV a. 

Simulations reveal the vertex finding efficiency is max- 
imum for total charged particle multiplicity of order 5 
and greater in the TPC. We studied the event counting 
efficiency of this analysis with a simple simulation based 
on events generated with the HIJING model [s^, and 
found the event counting efficiency is maximum for pro- 
duced charged particle multiplicities (in the range -0.5 
< rj < 0.5) exceeding 12. The vertex counting efficiency 
is of order 90% for multiplicities larger than 5, and falls 
abruptly to zero for smaller values. For this reason, the 
analysis presented in this work is limited to reference 
multiplicities in excess of 10 and 17 for Au + Au and 
Cu + Cu collisions where it is deemed minimally biased 
or unbiased. 

In order to eliminate track splitting we restricted our 
analysis to charged particle tracks producing more than 
20 hits within the TPC where 50% of these hits were 
included in the final fit of the track. 



III. NET CHARGE FLUCTUATION RESULTS 

We present, in Fig. [TJ measurements of the dynamical 

net charge fluctuations, ^dyn, as a function of collision 

centrality in Au + Au collisions at yfSNN = 19.6, 62.4, 
130, and 200 GeV, Cu + Cu coUisions at = 62.4 

and 200 GeV. 

In Fig. [U we see that the dynamical net charge fluc- 
tuations, in general, exhibit a monotonic dependence on 
the number of participating nuclcons. At a given num- 
ber of participants the measured fluctuations also ex- 
hibit a modest dependence on beam energy, with dyn 

magnitude being the largest in Au + Au collisions at 

y/sMN — 19.6 GeV. The dyn values measured for 

p + p collisions at ^/s — 200 GeV amounts to -0.230 ± 
0.019(stat). 
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FIG. 1: (Color online) Dynamical net charge fluctuations, 
!/+_,dyn, of particles produced within pseudorapidity jryj < 
0.5, as function of the number of participating nucleons. 



We first discuss the energy dependence of the fiuctu- 
ations. The collision centrality dependence is addressed 
in the following section. 



A. Beam Energy and Size Dependence 

A study of the net charge fluctuation dependence on 
the beam energy is of interest given that it can potentially 
reveal a change in the magnitude of the fluctuations and 
signal the formation of QGP. 

We conduct this study primarily on the basis of the 
0-5% and 0-10% most central collisions in Au -\- Au and 
Cu -f Cu collisions, respectively. Extensions to less cen- 
tral and peripheral collisions are possible but subject to 
additional uncertainties raised by small systematic errors 
involved the collision centrality determination. 

As already stated in the Introduction, charge conser- 
vation and the finite size of the colliding system intrinsi- 
cally limit the magnitude of the net charge correlations. 
Intuitively, one expects charge conservation effects to be- 
come progressively smaller with increasing charged parti- 
cle multiplicity. Charge conservation effects are nonethe- 
less definite at all beam energies and produced multiplici- 
ties. Specifically, one estimates that charge conservation 

implies a minimum value of order ^dyn = ^4/A'47r, 

where N^.,^ is the total charged particle multiplicity pro- 
duced over 47r (see for a derivation of this esti- 
mate). This estimate was obtained [1^ assuming that 
charge conservation implies global correlations but no 
dependence of these correlations on rapidity. Therefore, 
charge conservation effects may be different than those 
estimated in this work. Nonetheless, for simplicity, we 
use the above expression to estimate the effects of charge 
conservation on the dynamical net charge fiuctuations. 

Corrections to vj^ ^dyn for system size and charge con- 
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servation require knowledge of the total charged par- 
ticle multiplicity. Although, strictly speaking, no ex- 
periment at RHIC actually measures particle produc- 
tion with complete coverage, the PHOBOS experiment 
comes the closest with a rapidity coverage of \ri\ < 5.4 
over 2tt azimuthal angles and a minimum transverse mo- 
mentum of order 100 MeV/c. PHOBOS has published 
data on total measured charged particle multiplicities of 
Au + Au coUisions at _Js^ = 19.6, 62.4, 130 and 200 
GeV [13, m, m, EO, Ell, lH and Cu + Cu collisions at 
= 62.4 and 200 GeV We infer charged par- 
ticle multiplicities for p + p collisions at y/s = 200 GeV 
based on charged particle multiplicity per participant re- 
ported by PHOBOS We correct for differences in 
collision centralities between the PHOBOS and STAR 
measurements using a linear interpolation based on the 
two most central bins measured by PHOBOS. Number of 
participating nucleons (Npart), total multiplicities (Nch), 

uncorrected {v^ ^dyn) and corrected values (i^+°I^(jyn) of 

dyn are shown in Table [T] for p + p collisions at ^/s — 

200 GeV, all four energies in Au + Au collisions and two 
energies in Cu + Cu collisions. 



TABLE I: Number of participating nucleons, total multiplic- 
ity, uncorrected and corrected ,dyn values for p + p colli- 
sions at ^/s — 200 GeV, four energies in Au + Au collisions 
and two energies in Cu + Cu collisions. 



System & Energy 


part 






corr 


p + p 200 GeV 


2 


22 


-0.2301 


-0.04407 


Au + Au 200 GeV 


351 


5092 


-0.0024 


-0.00163 


Au + 130 GeV 


351 


4196 


-0.0021 


-0.00121 


Au + Au 62.4 GeV 


348 


2788 


-0.0029 


-0.00146 


Au + Au 19.6 GeV 


348 


1683 


-0.0035 


-0.00113 


Cu + Cu 200 GeV 


98 


1410 


-0.0071 


-0.00430 


Cu + Cu 62.4 GeV 


95 


790 


-0.0093 


-0.00437 



The corrected i^'^I7^yn values of the dynamical net 
charge fluctuations are shown in Fig. [2] as function of 
beam energy for 0-5% central Au + Au collisions with 
solid squares (in red color online) and for 0-10% central 
Cu + Cu collisions with solid circles (in black online). 
The displayed error bars include (a) the statistical errors 
involved in the measurement of ^dyn and (b) the to- 
tal charged particle multiplicities. The boxes show our 
estimates of the systematic errors involved in the mea- 
surements of both quantities. Data from this work are 
compared to corrected dynamical net charge fluctuations 
values by the PHENIX and CERES collaborations. The 
PHENIX point (triangle, in blue color online) is calcu- 
lated (as already discussed in JJ) from data published on 
the basis of the ljq observable Q and corrections based 
on total multiplicities measured by PHOBOS (as per val- 
ues shown in Table The CERES data points (star, in 
black online), obtained for Pb + Au collisions, are ex- 
tracted from their published results [1^. They include 
estimates of the systematic errors (open rectangles) as 
well as statistical errors (solid lines). 

We first note that the PHENIX and STAR points mea- 
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FIG. 2: (Color online) Corrected values of dynamical net 
charge fluctuations (i'^°r^tjyn) as a function of y/SNN- See 
text for details. 



sured at 130 GeV are in quantitative agreement as al- 
ready reported 0]. The large error bar associated with 
the PHENIX measurement stems mainly from systematic 
uncertainties associated with corrections for detection ef- 
ficiencies 0. We observe additionally that the STAR 
19.6 GeV measurement is in agreement with a measure- 
ment by CERES at the same energy. The STAR mea- 
surements in Cu + Cu collisions show a sharp increase 
in magnitude. This difference could partly be attributed 
to the difference in the number of participating nucleons 
in Au + Au and Cu + Cu collisions at 0-5% and 0-10% 
centralities, respectively. However, the magnitude of cor- 
rected dynamical fluctuations in Cu + Cu collisions when 
scaled by the ratio of number of participants in Cu + Cu 
collisions to number of participants in Au + Au collisions 
is - 0.0009 ± 2xl0"5(stat) ± 6xl0"5(sys) and -0.001 ± 
2xl0"5(stat) ± 8xl0-^(sys) at = 62.4 and 200 

GeV, respectively. We also note that CERES reports a 

dramatic reduction in the magnitude of ^dyn at the 

lowest energy measured at SPS. We thus conclude that 
net charge fluctuations corrected for charge conservation 
show no obvious beam energy dependence in the range 
from 19.6 to 200 GeV. However, there is a clear system 
size dependence when comparing Au + Au to Cu + Cu 
collisions. 

Below 19.6 GeV there appears to be a decrease in the 
magnitude of I'^iydyn the lowest SPS energies. Dif- 
ference between STAR and CERES results may in part 
stem from differences in pseudorapidity acceptance. 

Measurements at the SPS have shown that particle 
production at 5 GeV and lower energies is dominated by 
baryons while meson and resonance production become 
increasingly dominant at energies above 19.6 GeV. This 
suggests that the change in dynamical net charge fluc- 
tuations below 19.6 GeV might, in part, be due to this 
shift in particle production dominance. It is also con- 
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ceivable that the differences between the values measured 
below and above 19.6 GeV may result from changes in 
the collision dynamics and final state interaction effects 

[H E [S Q m, K [13, E US mi • 



B. Collision Centrality Dependence 

The observed monotonic reduction of the magnitude 

of _dyn with increasing number of participants, seen 

in Fig. [U arises principally from the progressive dilution 
of two- particle correlation when the number of particle 

sources is increased. In fact, one expects ^dyn to be 

strictly inversely proportional to the number of partici- 
pating nucleons or the produced particle multiplicity if 
Au -\- Au collisions actually involve mutually indepen- 
dent nucleon-nucleon interactions, and rescattering ef- 
fects may be neglected. 

We investigate the possibility of such a scenario by 
plotting the dynamical fluctuations scaled by the mea- 
sured particle multiplicity density in pseudorapidity 
space {dNch/di]) in Fig. [3fa). Data from Au + Au col- 
lisions at various energies are shown with solid symbols 
while data from Cu + Cu collisions at 62.4 and 200 GeV 
are shown with open symbols. Values of dNch/drj used 
for the scaling correspond to efficiency corrected charged 
particle multiplicities measured by STAR f46'| and PHO- 
BOS [37, 38, 39, ,40, 41,i2,i43|]- We note that the correc- 
tion applied in Section A to account for charge conserva- 
tion, is useful to study the energy dependence of the net 
charge fluctuations. Its use for centrality, pseudorapidity, 
and azimuthal dependencies is, however, not warranted 
given that insufficient data are available to reliably ac- 
count for charge conservation effects. Also, the applied 
correction is model dependent, i.e., assumes charge con- 
servation applies only globally [29| . 

We note from Fig. [3]Ja) that the magnitude of ly^ ,dyn 

scaled by dNch / drj for Au -I- Au 200 GeV data is differ- 
ent from the rest of the data. This could partly be at- 
tributed to the larger multiplicity produced in Au + Au 
200 GeV. We additionally observe that all four distri- 
butions exhibit the same qualitative behavior: the am- 
plitude ,dyndNch/dr]\ is smallest for peripheral col- 
lisions, and increases monotonically by ^^40% in central 
collisions in Au + Au and Cu + Cu systems. The ob- 
served ,dyndN ch / dri\ increases with the increase in 

collision centrality. The dashed line in the figure corre- 
sponds to charge conservation effect and the solid line 
to the prediction for a resonance gas. The figure in- 
dicates that dynamical net charge fluctuations, scaled 
by dNch /drj are rather large. Most central collisions in 
Au+Au 200 GeV approach the prediction for a resonance 
gas Indeed, observed values of vj^ ^dyn are inconsis- 
tent with thosepredicted based on hadronization model 
of Koch et al. [J 0, [1] . Given recent observations of el- 
liptic flow, suppression of particle production at high px 
(Raa 0-2), and two-particle correlation functions in- 
dicating the formation of a strongly interacting medium 



(sQGP) in A + A collisions at RHIC energies, this sug- 
gests that the signal predicted by the authors [J [3, Q 
may be washed out by final state interactions, diffusion, 
expansion, collision dynamics, st ring fusion [47] or other 
effects [iil, [H, [H, H [H [il M, mhh., M, M , some of 
which were discussed in the introduction. 

Changes in the collision dynamics with increasing cen- 
trality are indicated by these data. Such a conclusion 
should perhaps not come as a surprise in view of the 
large elliptical flow, and the significant reduction of par- 
ticle production at high transverse momenta reported 
by aU RHIC experiments We also note the PHO- 

BOS collaboration has reported that the charged par- 
ticle multiplicity per participant nucleon pair rises sub- 
stantially with increasing number of participants. They 
report a value of dNch/ drj/ {N part/ "2) of order 3.9 in cen- 
tral 200 GeV Au -I- Au collisions compared to a value 
of 2.5 in p -\- p collisions at the same energy ^39!]. This 
amounts to a 56% increase, similar in magnitude to 

that of ,dyndNch/drj\ measured in this work. We 

thus infer that much of the centrality dependence of 
W+-,dyndNch/dri\ is due to the rise of dNch/ drj/ {Npart/"^) 
with increasing Npart- 

In order to validate this assertion, we plot in Fig. [S^b) 
the dynamical fluctuation scaled by the number of par- 
ticipants, Nparti^-\ ,dyn as a function of the number of 

participants. Vertical error bars represent statistical un- 
certainties. Values of NpartV-\ ,dyn exhibit a small de- 
pendence on the collision centrality at all four measured 
energies in Au+Au collisions and two energies in Cu+Cu 
collisions. The measured data scaled by the number of 
participants {Npart) are thus consistent with either no or 
a very weak centrality dependence. However, a definite 
system size and energy dependence is observed. This im- 
plies that the strength of the (integrated) net charge two- 
particle correlation per participant exhibits essentially no 

dependence on collision centrality. We also scale ^dyn 

with the number of binary collisions, shown in Fig. OJc). 
While we observe that the datasets follow a common 
trend, dyn clearly exhibits dramatic collision central- 
ity dependence. Such a dependence is, however, expected 
given that the measured dynamical net charge fluctua- 
tions are dominated by low momentum particles with 
large cross-section for which binary scaling does not ap- 
ply. The statistical errors on dyn and the scaling 

factors used in Fig. [31^ a), [3{b) and are added in 
quadrature. 



C. Longitudinal and Azimuthal Dependencies of 
the Dynamical Fluctuations 

Pratt et al. [i^, HO] have argued that the width of 
longitudinal charge balance functions should significantly 
narrow in central Au -\- Au collision relative to peripheral 
collisions or p -f p collisions due to delayed hadroniza- 
tion following the formation of a QGP. STAR has in fact 
reported that, as predicted, a narrowing of the balance 
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function does occur in central Au + Au collisions rel- 
ative to peripheral collisions [lH. We note, however, 
as already pointed out by Pratt et al. and more re- 
cently by Voloshin p^, radial flow produced in heavy 
ion collisions induces large position-momentum correla- 
tions which manifest themselves in angular, transverse 
momentum, and longitudinal two-particle correlations. 
The observed narrowing of the longitudinal charge bal- 
ance function therefore cannot be solely ascribed to de- 
layed hadronization. It is thus important to gauge the 
change in two-particle correlations imparted by radial 
flow effects. As a first step towards this goal, we present 
studies of the net charge fluctuation dependence on the 
integrated pseudorapidity and azimuthal ranges. 
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FIG. 3: (Color online) Dynamical net charge fluctuations, FIG. 4: (Color online) Dynamical fluctuations !^+-,dyn, nor- 
!/+_,dyn, of particles produced with pseudorapidity \ri\ < 0.5 malized to their value for |t?| < 1, as function of the inte- 
scaled by (a) the multiplicity, dNch/drj. The dashed line cor- grat ed ps eudorapidity range, (a) data for Au + Au coUisions 
responds to charge conservation effect and the solid line to the at ^snn = 62 4, 20 GeV (0-5%) along with data for Cu+Cu 
prediction for a resonance gas, (b) the number of participants, coUisions at Vsatjv = 62.4, 200 GeV (0-10%), are compared 
and (c) the number of binary coUisions. to inclusive p + p data at ^ = 200 GeV, and (b) data for 

Au + Au coUisions at = 62.4, 200 GeV (30-40%) along 

with data for Cu + Cu collisions at y^sjvjv = 62.4, 200 GeV (0- 
10%), are compared to inclusive p + p collision data at ^/s — 
200 GeV. 
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We plot in Fig. \^a,) values of ,dyn(T?) measured for 

different ranges of pseudorapidity, rj. In order to com- 
pare data measured at different centralities, beam ener- 
gies and system size, measured values are normalized by 

the magnitude of i/^ ,dyn(?7) for a pseudorapidity range 

I77I < 1 {v-i ,dyn(l))- The data shown in Fig. ID^a) are 

from Au + Au collisions at y^sjvw = 62.4 and 200 GeV, 
Cu + Cu collisions at 62.4 and 200 GeV, and p-^ p data 
obtained at 200 GeV. One finds the magnitude of the 
normalized correlation is maximum for the smallest pseu- 
dorapidity ranges and decreases monotonically to unity, 
at all energies and centralities, with increasing pseudora- 
pidity range. 

The dynamical fluctuations being essentially a mea- 
sure of two-particle correlation dominated by the 

term, one finds, as expected, that the correlation is 
strongest for small rapidity intervals, and is increasingly 
diluted (reduced) for larger intervals. For example, in 
Cu + Cu collisions at ^/snn — 200 GeV the typical val- 
ues of and are 0.99256, 0.992518 and 
0.996099, respectively. One observes that the magni- 
tudes of ,dyn(?7)/l^H ,dyn(l)| m Cu -\- Cu COllisionS 

at 62.4 and 200 GeV are quite different from Au + Au 
collisions at comparable energies. This shows that the 
collision dynamics in p + p collisions, 0-10% Cu -\- Cu 
and 0-5% Au -\- Au collisions are significantly different. 
Indeed, we find the relative magnitude of the correla- 
tions measured for |?7| < 0.5 increases by nearly 25% for 
Au + Au 200 GeV relative to those mp-\-p. Note in par- 
ticular that the slope {dv^ ,dyn/dr]) inp+p, Cu+Cu and 

Au + Au systems depends on the correlation length (in 
pseudorapidity): the shorter the correlation, the larger 
the slope. The observed distributions then indicate that 
the correlation length is shorter for central collisions and 
for larger systems, in agreement with the observed re- 
duction of the charge balance function [5^. The larger 
values of the slopes observed for most central collisions 
(as well as for larger systems) indicate correlated pairs 
of negative/positive particles tend to be emitted closer 
in rapidity than those produced in peripheral Au + Au 
OT p + p collisions. Authors of Ref. 49] have proposed 
that a reduction of the width of the balance function, and 
conversely a relative increase of short range (jryj < 0.5) 
correlations, could signal delayed hadronization. The ob- 
served increase in the correlation, reported here, might 
however also result from the strong radial flow believed 
to exist in central Au + Au collisions. 

A comparison of Au + Au collisions at ^/snn — 62.4, 
200 GeV (30-40% central) is made with Cu + Cu colli- 
sions at the two energies for 0-10% centrality in Fig. Hfb), 
as these centralities correspond to approximately same 
number of participant nucleons. We observe that the 
magnitude of normalized correlation is similar for both 
systems at the same beam energy, thereby suggesting 
that the magnitude and the width of the charge particle 
correlation depends mainly on the number of participants 
and collision energy but little on the colliding systems. 

To understand the role of radial flow in net charge fluc- 



tuations measured in a limited azimuthal range (i.e. less 

than 27r), first consider that the magnitude of ,dyn 

is, in large part, determined by the abundance of neu- 
tral resonances (such as the p{77Q)). The decay of neu- 
tral resonances into pairs of charged particles increases 
the charged particle multiplicity without affecting the 
variance of the net charge. An increasing fraction of 
neutral resonances (relative to other particle production 
mechanisms) therefore leads to reduced magnitude of 
^'-1 ^dyn- Consider additionally that large radial flow ve- 
locity should lead to a kinematical focusing of the decay 
products in a narrow cone. The opening angle of the cone 
will decrease with increasing radial velocity boost. One 
thus expects that while measuring i/j, ^dyn m a small az- 
imuthal wedge, one should have greater sensitivity to the 
level of kinematical focusing, i.e. the magnitude of the 
dynamical net charge fluctuation (correlation) should in- 
crease with the magnitude of the radial flow velocity. Az- 
imuthal net charge correlations should therefore be rather 
sensitive to the magnitude of the radial flow velocity. 

Fig. [SJa) and (b) display azimuthal net charge corre- 
lations integrated over azimuthal angle ranges from 10 
to 360 degrees for Au + Au and Cu + Cu collisions at 
200 GeV. An azimuthal wedge of, for example, 90 de- 
grees would divide the complete phase space into four 
sectors, where we denote each sector as a bin. The fig- 
ure shows results from nine azimuthal wedges obtained 
after averaging z^+.^dyn values for all bins in each wedge. 
The errors shown in EJa) and (b) show the statistical 
errors of the averaged values for each wedge size. We 
also verified that for small wedge angles (e.g. 90 degrees 
and smaller), the variances of the measured values, for 
wedges of a given size, have a magnitude similar to the 
errors of the averages. Data are shown for seven colli- 
sion centrality bins in Au + Au collisions in EI a) and for 
five centrality bins in Cu + Cu collisions in[5^b). Note 
that the absolute magnitude of the correlation decreases 
from the most peripheral to the central collisions as a 
result of progressive dilution with increasing number of 
participants. The variation of the shape of the correla- 
tion function with the size of the azimuthal acceptance 
is of greater interest. One finds the correlation functions 
measured in the most central collisions decrease mono- 
tonically in magnitude with increasing azimuthal wedge 
size whereas they exhibit a more complicated behavior 

for most peripheral collisions. One expects vj^ ^dyn to be 

rather small for very small acceptance (azimuthal wedge) , 
i.e., when the size of the acceptance is smaller than the 

typical correlation length. This explains why dyn| 

decreases sharply for small angles in peripheral collisions. 
It is remarkable, however, to note that this behavior is 
not observed in most central collisions with the angular 
ranges considered thereby indicating a change in the par- 
ticle correlation length qualitatively consistent with the 
reduction of the balance function in central collision 

aheady reported by STAR 

Fig. [5] shows a comparison of Au+Au and Cu+Cu col- 
lisions at similar number of participating nucleons. The 
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FIG. 5: (Color online) Dynamical fluctuations ,dyn, as 

a function of the integrated azimuthal range (f) for selected 
collision centralities for (a) Au + Au collisions at ^/smn ~ 
200 GeV, and (b) Cu + Cu collisions at y/s^ = 200 GeV. 
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FIG. 6: (Color online) Dynamical fluctuations ,dyn, as 

a function of the integrated azimuthal range 4> for similar 
number of participating nucleons for Au + Au and Cu + Cu 
collisions at s/snn ~ 200 GeV. 



magnitude of ^dyn with respect to azimuthal angle, 

(/>, is similar for similar number of participating nucleons 
in both systems with best agreement for collisions with 
more than 20 participants. The agreement for the most 
peripheral collisions studied is weaker, but we speculate 
that vertex inefficiencies and fluctuations in the number 
of participants should account for this weaker agreement. 
We also observe a change in shape with centrality. Both 
systems show similar monotonic dependence for central 

collisions, whereas, the magnitude of v^^ .dyn reaches a 

maximum for a small azimuthal wedge for peripheral col- 
lisions. The error bars shown here are statistical only. 



IV. SYSTEMATIC ERROR STUDIES 

While ^dyn is a robust observable and shown to 

exhibit essentially no dependence on efficiencies, it may 
nonetheless be subject to limited systematic effects as- 
sociated with the measurement process. We investigated 
dependencies on the longitudinal position of interaction 
vertex (z- vertex), the effect of resonance feed downs, 
event pile-up, track reconstruction and pt resolution. 

The dependence of ^dyn on the longitudinal posi- 
tion of the interaction vertex might arise because of the 
restricted acceptance of the TPC on which these analyses 

are based. We thus measured ,dyn by binning events 

according to the z-vertex in steps of 5 cm for positions 
varying in ranges 5 < \z\ < 30 cm and found deviations 
in dyn to be 1% or less. 

The dyn measurement presented in this paper is 

meant to be representative of particles produced by 
Au + Au, Cu + Cu or p+p collisions. By design, one thus 
seeks to eliminate effects from secondary decays (e.g., 
K p + 7r~) or secondary particle production within 
the detector. This is accomplished by limiting the anal- 
ysis to tracks that appear to originate from the collision 
vertex. Indeed a cut of track distance-of-closest approach 
(DCA) to the collision vertex with a value of 3 cm is used 
to select primary particles and reduce those produced by 
decays and secondary interactions. The large value of 
DCA used in this analysis is due to finite DCA resolu- 
tion and is also intended to maintain large track detec- 
tion efficiency, which is needed especially for the ,dyn 

analysis with respect to the longitudinal and azimuthal 
acceptance. However, with a large value of the DCA 
cut, one ends up counting particles produced by weak- 
decays (e.g., A or K'^) as primary particles. In particular, 
with kaons {K^) representing a small fraction of all neu- 
tral particles produced, one expects pions from decays of 
these particles to increase the accepted charged particle 
multiplicity but with only a minor impact on the vari- 
ance of the measured net charge. This implies ,dyn 

should be subject to a systematic decrease in magnitude 
when accepting weak-decay feed down. We thus studied 

_dyn for smaller DCA cuts of 2 cm and found _dyn| 

decreases by roughly 1% at all collision centralities. 
and A have a decay length in excess of 2.7 cm. Given the 
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rather limited resolution of the measurement, the DC A 
of the decay products is spread to values over a range 
extending more than 3 cm and, thereby, form a modest 
background to the primary particles. Assuming the con- 
tributions of Kg and A are roughly uniform within the 
3 cm DCA cut considered, we expect that a 2 cm DCA 
cut reduces the background by approximately 30%. We 
observe this change of cut leads to a 1% reduction in the 

magnitude of z/^ dyn- We thus conclude that and A 

contamination amounts to a contribution of a few percent 
only. 

Another important source of secondary tracks not 
completely eliminated by the DCA cut arc elec- 
trons/positrons. While a finite electron primary yield is 
expected from decays of D-mcsons and B-mcsons, from 
Dalitz decays of 7r° and rj, the bulk of electrons/positrons 
observed in the TPC are from secondary interactions 
leading to pair production, and Compton photo-electron 
production. Elimination of electrons/positrons is, in 
principle, partly achievable based on cuts on track 
dE/dx. However, because electrons and pions of low mo- 
menta experience similar energy loss in the TPC gas, a 
cut on the track dE/dx also eliminates a large amount 
of pions thereby effectively creating a "hole" in the pion 
acceptance (with respect to their momentum). We thus 
carried out the analysis reported in this paper by in- 
cluding the electrons/positrons. Again in this case, since 
electrons and positrons are typically created in pairs, this 
may lead to an increase in the integrated charged parti- 
cle multiplicity with little impact on the net charge vari- 
ance. One thus expects inclusion of the electrons should 

produce a systematic shift in the magnitude of ^dyn- 

To verify this we carried out a measurement of _dyn 

when electrons (and consequently also pions) are elim- 
inated on the basis of dE/dx cut. The dE/dx cut is 
accomplished using the truncated mean of the measured 
dE/dx samples along the track and the track momentum. 
Tracks were excluded whenever the measured dE/dx fell 
within two standard deviations of the mean value ex- 
pected for electrons of a given momentum. We found 

that when electrons are eliminated, _dyn| increases 

by as much as 3.5% in magnitude. This shift may how- 
ever not be entirely due to the suppression of electrons. 
Indeed, by cutting electrons, one also reduces pion ac- 
ceptance in transverse momentum. We have reported in 
Section III.C. that v+^^^yn exhibit a modest dependence 
on the size of integrated longitudinal and azimuthal ac- 
ceptances. However, a similar (but weaker) dependence 
on the transverse momentum is expected. It is thus plau- 
sible the shift by 3.5% may in part result from a reduction 
of pion acceptance. Electron contamination is thus con- 
sidered a source of systematic error of the order of 3.5% 
in our measurement of ,dyn- 

Au -f- Au and Cti + Cu data acquired during nms IV 
and V were subject to pile- up effects associated with 
large machine luminosity obtained during those years. 
The pile-up may result in two collisions being mistaken 
as one and treated as such, thereby leading to artificially 



large multiplicities and increased variances. Therefore, in 
order to reject pile-up events, dip angle cuts were intro- 
duced in the present analysis. The dip angle is defined as 
the angle between the particle momentum and the drift 
direction, 9 = cos~^ (pz/p). The dip angle cut is based on 
the average dip angle distribution of all tracks in a given 
event. We found the dip angle is correlated with the 
vertex position and features a width distribution which 
is Gaussian at low luminosities. We thus reject pile-up 
events that are beyond two standard deviations of the 
mean of the distribution for a particular centrality and 

vertex position. We found y^ ^dyn changes by less than 

1% when the dip angle cut is used. 

We also checked the effect of efficiency variation within 
the acceptance of interest. The efficiency is known in par- 
ticular to progressively reduce from a maximum value for 
Pt > 200 MeV/c to zero for px <100 MeV/c. We deter- 
mined an upper bound of the effect of pr dependence by 

measuring i^^ ^dyn with px thresholds of 150 MeV/c and 

200 MeV/c. We found changes of dyn are typically 

negligible within the statistical accuracy of our measure- 
ment and amount to at the most 1.5%. 

Total systematic error contribution increases from 8% 
to 9% from central to peripheral collisions in Au + Au 
and Cu + Cu collisions at ^/snn = 200 GeV. Similarly, 
systematic errors amount to 8% in peripheral collisions 
and 7% in central collisions in Au + Au and Cu + Cu 
collisions at ^/sNN — 62.4 GeV. The systematic errors 

on ,dyn from different sources mentioned above are 

added linearly. 



V. SUMMARY AND CONCLUSIONS 

We have presented measurements of dynamical net 
charge fluctuations in Au + Au collisions at y^SNN = 
19.6, 62.4, 130, 200 GeV, Cu + Cu collisions at ^/s]^ = 
62 A, 200 GeV and p + p collisions at ^ = 200 GeV, 
using the measure y^ ^dyn- We observed that the dy- 
namical net charge fluctuations are non vanishing at all 
energies and exhibit a modest dependence on beam en- 
ergy in the range 19.6 < ^/sjvjv < 200 GeV for Au+Au as 
well as Cu + Cu collisions. Dynamical fluctuations mea- 
sured in this work are in quantitative agreement with 
measurements by the CERES collaboration at ^/sNN = 
17.2 GeV and PHENIX collaboration at ^/s]^ = 130 
GeV. However, measurements by CERES at lower beam 
energy (<17.2 GcV) exhibit much smaller dynamical net 
charge fluctuations perhaps owing to a transition from 
baryon to meson dominance in the SPS energy regime. 
We also found the dynamical net charge fluctuations vi- 
olate the trivial 1/Nch scaling expected for nuclear col- 
lisions consisting of independent nucleon-nucleon inter- 
actions. However, one finds that dyn scaled by the 

number of participants exhibits little dependence on col- 
lision centrality but shows modest dependence on colli- 
sion systems. Measured values of ^dyn are inconsis- 
tent for all systems and energies with the predictions of 
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the QGP hadronization model of Koch et al. P, 0, Q. 
Given the reported observations of a strongly interact- 
ing medium m A + A collisions at RHIC, this suggests 
that the assumptions of the hadronization by Koch et 
al. are invalid, or that some final state interaction pro- 
cess washes out the predicted signal. Scaled dynamical 

net charge fluctuations ,dyndN ch I dri\ grow by up to 

40% from peripheral to central collisions. We speculated 
that the centrality dependence arises, in part due to the 
large radial collective flow produced in Au + Au collisions 
and proceeded to study fluctuations as a function of az- 
imuthal angle and pseudorapidity. Our analysis showed 
dynamical fluctuations exhibit a strong dependence on 
rapidity and azimuthal angular ranges which could be 
attributed in part to radial flow effects. 
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